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Abstract 

We present a detailed analysis of the role of the magnetic dipole-dipole interaction in cold and 
ultracold collisions. We focus on collisions between magnetically trapped NH molecules, but the 
theory is general for any two paramagnetic species for which the electronic spin and its space-fixed 
projection are (approximately) good quantum numbers. It is shown that dipolar spin relaxation is 
directly associated with magnetic-dipole induced avoided crossings that occur between different adi- 
abatic potential curves. For a given collision energy and magnetic field strength, the cross-section 
contributions from different scattering channels depend strongly on whether or not the corre- 
sponding avoided crossings are energetically accessible. We find that the crossings become lower in 
energy as the magnetic field decreases, so that higher partial-wave scattering becomes increasingly 
important below a certain magnetic field strength. In addition, we derive analytical cross-section 
expressions for dipolar spin relaxation based on the Born approximation and distorted-wave Born 
approximation. The validity regions of these analytical expressions are determined by comparison 
with the NH + NH cross sections obtained from full coupled-channel calculations. We find that the 
Born approximation is accurate over a wide range of energies and field strengths, but breaks down 
at high energies and high magnetic fields. The analytical distorted-wave Born approximation gives 
more accurate results in the case of s-wave scattering, but shows some significant discrepancies for 
the higher partial- wave channels. We thus conclude that the Born approximation gives generally 
more meaningful results than the distorted-wave Born approximation at the collision energies and 
fields considered in this work. 
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I. INTRODUCTION 



The ability to produce and trap atomic and molecular species at sub-kelvin temperatures 



offers numerous exciting possibilities in condensed-matter physics [l 
7|, high-precision spectroscopy js iJl, and physical chemistry |l3 



uantum computing 
22[. Since the experi- 



23|, |2J], major advances have been 



mental realization of the first Bose-Einstein condensates 
made in the field of ultracold atomic gases. It is now well established that alkali-metal atoms 
can be efficiently cooled into the ultracold regime using a combination of laser cooling and 
evaporative cooling. However, laser cooling is not applicable to all atomic species, and is 
particularly difficult for molecules {25I. In the last few years, several methods have been 
developed that aim at producing (ultra^cold molecular gases at relatively high densities. 
Techniques such as photoassociation [ 



and magnetic Feshbach association [27| employ 
an indirect scheme in which ultracold molecules are formed by pairing up pre-cooled, ul- 
tracold atoms. These methods are, however, currently limited to molecules consisting of 
;wo alkali-metal atoms. Direct-cooling methods such as Stark 28|] and Zeeman deceleration 
29 1 . molecular-beam guiding 30 1, and buffer-gas cooling [si^ apply to a much wider range 
of molecular species, but require a second-stage cooling technique to reach the ultracold 
regime. Although several theoretical studies have shown that e.g. sympathetic cooling of 



cold molecules with ultracold co-trapped atoms 



37 



32|-|36| or molecular evaporative cooling 



is likely to be successful, this is yet to be demonstrated experimentally. 
Second-stage cooling methods such as forced evaporative cooling require strong elastic 
collisions that thermalize the gas cloud as the trap depth is slowly reduced 40|. Inelastic 
collisions, in which the internal quantum state of at least one of the collision partners is 
changed, can induce heating of the gas and trap loss. A detailed understanding of the 
interparticle interactions that govern these inelastic processes is thus crucial for assessing 
the feasibility of second-stage cooling. One of the most important inelastic loss mechanisms 
for trapped paramagnetic species is dipolar spin relaxation, which arises from the magnetic 
dipole-dipole interaction between the magnetic moments of the particles. For many spin- 



polarized atomic gases such as hydrogen 40 1, lithium j41|, nitrogen 42 1, and chromium 



43|, 



3ut also for atom-mo^ 
M, and NH + NH 



ecu 



38 



e and molecule-molecule systems such as Li + NH 3J|, N -|- NH 



391], the interparticle dipolar spin-spin interaction is indeed the 



dominant source of trap loss. 
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In this paper, we provide a comprehensive study on the role of the magnetic dipolar 
interaction in cold and ultracold collisions. Specifically, we consider collisions between mag- 
netically trapped bosonic ^^NH(X^S~) molecules, but the theory should be general for any 
(ultra) cold paramagnetic species. We assume that the molecules are in their vibrational 
and rotational ground states, as is the case experimentally 16|. For NH + NH, there are 
three spin-changing mechanisms that can induce trap loss: the intramolecular spin-spin and 



spin-rotation couplings, and the intermolecular magnetic dipolar coupling term j45|. Pre- 
vious theoretical work [sqI] has shown that the intermolecular magnetic dipole interaction 
is the main spin-relaxation mechanism for NH-NH at low collision energies and small to 
moderate magnetic field strengths. It was also shown, in the same paper, that the dipolar 
spin-spin coupling term induces certain avoided crossings between different adiabatic poten- 
tial curves, which in turn give rise to spin-changing transitions. That is, the spin-fiip due 
to the intermolecular magnetic dipolar interaction can be qualitatively understood in terms 
of the avoided curve crossings 39|]. In the present work, we discuss the infiuence of these 
crossings on the cross section in much greater detail. We also provide analytical expressions 
for the dipolar spin-relaxation cross section based on the Born approximation (BA) and 
distorted- wave Born approximation (DWBA). We compare the analytical results with the 
cross sections obtained from rigorous close-coupling (CC) calculations, and show that the 
results are in excellent agreement over a wide range of collision energies and magnetic field 
strengths. 

This paper is organized as follows. In Sec. Ill Al we briefiy describe the details of the 
CC calculations. The derivations of the BA and DWBA cross sections are given in Sees. 
IIIBI and III CI respectively, and the results are discussed in Sec. IIIII The numerical results 
are presented in Sec. IIII At with a particular emphasis on the role of the avoided curve 
crossings, and the validity of the analytical BA and DWBA cross sections is detailed in Sec. 
IIIIBI Finally, concluding remarks are given in Sec. IIVI 



II. THEORY 

Throughout this paper, we will focus on collisions between two bosonic ^^NH(X^S~) 
molecules in their magnetically trappable, low-field seeking states \Sa = 1,^5a ~ ^)\Sb = 
1, Msg = 1). Here Si denotes the total electronic spin of the monomers {i = A, B) and M5. 
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is the spin projection onto the magnetic field axis. A colhsion complex of two such molecules 
is in the high-spin quintet 15 = 2, Ms = 2) state, with S denoting the total spin and Ms its 
space-fixed projection. Collisions that change either the Ms quantum number of the quintet 
state or the total spin S to yield singlet {S = 0) or triplet {S = 1) complexes will lead to 
immediate trap loss. 



A. Coupled-channel calculations 

In order to obtain numerical values for the collision cross sections of NH + NH, we have 
performed full CC calculations as a function of energy and magnetic field. The details of 



these calculations are given elsewhere 



39| and we provide only a brief description here. The 



NH-NH scattering Hamiltonian is written as 

fi^ d"^ P' .... 

where /i is the reduced mass of the complex, R is the intermolecular vector that connects the 
centers of mass of the monomers, R = \R\, L"^ is the angular momentum operator associated 
with rotation of R, V{R, coa, i^b) is the potential-energy surface for the quintet {S = 2) state 
of NH-NH, ua and ub describe the orientation of monomers A and B, Vmagn.dip(-R, Sa, Sb) 
is the intermolecular magnetic dipolar interaction between the two spins, and Ha and Hb 
are the Hamiltonians of the individual monomers. The magnetic dipole-dipole term is given 
by 

K„agn.dip(i«, Sa, Sb) = -VQqWb^ 5^(-l)^C2,_,(fi«)[5A ® 5^]^ , (2) 

where gs ~ 2.0023 is the electron (yf-factor, /iB is the Bohr magneton, a is the fine-structure 
constant, C2-q is a Racah-normalized spherical harmonic, Vtn = {Qr, ^r) describes the 
orientation of R, and the factor in square brackets is the tensor product of the monomer 
spin operators Sa and Sb- The monomer operators Hi correspond to the asymptotic molec- 
ular states and account for the monomer rotation, intramolecular spin-spin coupling, spin- 
rotation coupling, and Zeeman interaction. Hyperfine coupling is neglected. 

The scattering calculations were carried out in a symmetry-adapted basis set that ac- 
counts for the identical-particle symmetry of the system, 

1 

[2(1 + 5;,,J]V 



Z.sLaJ = 1^^7T-1 ^,h^^B)+v{-l)'\lBlA)]\LM,). (3) 
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Here r] defines the symmetry of tlie wave function with respect to molecular interchange, 
which is +1 for the bosonic ^^NH - ^^NH complex, e = (^—I^'^a+Nb+l jg ^j^g parity symmetry, 
which must be +1 for identical bosons in the same quantum state, and [7^, 7b) denotes the 



molecular rotation and spin functions in the space-fixed frame 



45| 



|7a7b) = \NaMn^)\SaMs^)\NbMn^)\SbMs,). (4) 

The basis set was truncated at Na = Nb = 2 and L = 6. Although this basis set is not 
fully converged, we have verified that the calculated cross sections are very similar to those 
obtained with A^^ = Nb = 3 in the region where the intermolecular dipole-dipole coupling 
is dominant, i.e. at ultralow energies and small to moderate field strengths. Increasing the 
rotational basis set does yield a larger cross-section contribution from the intramolecular 
spin-spin coupling, but this term becomes important only at energies above ~1 mK and 
fields above ~100 G. For a more general discussion on the issue of basis-set convergence, the 



reader is referred to Refs. 



38| and |39|. 



Let us now consider the identical-particle symmetry of the complex. Even though hyper- 
fine coupling is neglected, the symmetry of the nuclear-spin wave function should be taken 
into account when evaluating the exchange symmetry of the total wave function. We have 
assumed that both monomers are in their nuclear-spin stretched states (/ = Mj = 1), so 
that the nuclear-spin function is symmetric under exchange. Thus, we have 77 = +1 and 
e = +1. We also point out that, due to parity conservation, collisions between rotational 
ground-state molecules can only occur for even values of L. Furthermore, the conservation of 
the total angular momentum projection Ai = Mn^ + M^b + + + requires that 
any change in Mg^^ or Msg must be accompanied by a change in M^. It therefore follows 
that, in the ultracold regime, the s-wave spin-inelastic collision channel for magnetically 
trapped, rotational ground-state NH is dominated by the L = 2 outgoing partial wave. 

We performed the scattering calculations for each value of Ai and accumulated the re- 
sulting scattering ^-matrices to extract the cross sections. The calculations were carried 
out using a modified version of the MOLSCAT package 0, The propagation was per- 
formed using the hybrid log-derivative method of Alexander and Manolopoulos 481]. Prior 
to matching to asymptotic boundary conditions, an additional transformation was required 
to obtain the exact channel eigenfunctions j45|. This is because the intramolecular spin-spin 
coupling mixes states with Ni and Ni±2, which makes Ni, Mjsf., and M5- only approximately 



good quantum numbers. The exact molecular eigenstates will be denoted as 



17^7^) = \{Na,Sa)Ja,Mj^)\{Nb,Sb)Jb,Mj^). (5) 

We emphasize that the intramolecular coupling is relatively weak and Ni, M^^, and M5. 
may be treated as almost exact. Specifically, for the rotational ground state of ^^NH, the 
magnetically trapped component with Jj = Mj- = 1 contains 99.992% of |A^j = 0,Miy- = 
0,Si = l,Ms, = l). 



B. Born approximation 

In this section, we derive an analytical expression for the inelastic spin-changing cross 
section due to Knagn.dip based on the first-order Born approximation. This approximation 
assumes that the interaction between projectile and target is so weak that the initial and final 
states can be described by undistorted plane waves. We note that the BA has been previously 



used in the study of cold collisions in e.g. Refs. 



49|-|52| . The aim of the present work is to give 



a cross-section expression in closed form, and we therefore outline the complete derivation 
for the sake of clarity. The derived expression is general for any two paramagnetic species for 
which the electronic spin and its space-fixed projection are (approximately) good quantum 
numbers, e.g. for Hund's case (b) molecules and S'-state atoms, but we will apply it only to 
the case of NH(3S-) + NH(3S-). 

We start with the exact expression for the differential cross section (see e.g. Eq. (XIX. 19) 



of Ref. 



53|) 



dag^bi^a) _ 27r 1^,^ _(,,) |^(+) _(a)\|2 



f^Kfe.,7?\7l?1^nt|feW,7j\7lr)>|V(i?), (6) 



where a and b label the initial and final states, respectively, Qa = {(^kai'i^ka) ^"^^ = 
{6k,,,4'kt) describe the directions of the incoming and outgoing collision fiuxes, fci'*''' is the 
exact incident wave function with wavenumber ka, kb is a plane wave with wavenumber kb, 
7^^ and T"^"* denote the internal quantum numbers of the monomers for the initial and final 
states {i = a,b), V^nt is the interaction between the scattering particles, for which we take 
Vint = Knagn.dip, Va = fika/ is the velocity of the incident beam, and pb{E) = fikb/[h^{27iy] is 
the density of final states at energy E = ^fJ^v"^. The first-order Born approximation amounts 
to approximating the incident wave function as a plane wave, i.e. \ki'^^) ^\K)- The plane 
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waves are normalized to unit density and are mutually orthogonal, 

|fc) = e'*=-^, (7) 
{k\k') = {27vf6{k-k'). (8) 

Here S{k — k') represents the three-dimensional Dirac delta function. 

In the case of NH + NH, the asymptotic states "Ja and 7^ should be described as in 
Eq. ([5]). However, since we focus on collisions between rotational ground-state molecules, 
we may treat Ms^ and Msg as almost exact quantum numbers. Furthermore, taking into 
account that V^agn.dip acts only on the vector R and the electron-spin coordinates, we can 
omit the molecular rotational quantum numbers and write 

|7«> - |S^M«>, 

|7g^> - \SbMII). (9) 

The energies of the initial and final molecular states are now determined only by their 
Zeeman shifts. If we define the Zeeman levels relative to the initial state, the wavenumbers 
are ka = y/2jj,E/fi and ki, = ^J2^{E + qs^b^MsB) /Ti, where AMs is the total spin-change, 
AM5 = M^fj^ + M^g^ — Mgj^ — M^^, the term qsHb^MsB is the corresponding change in 
Zeeman energy, and B is the magnetic field strength. In the remainder of this paper, we 
will use B exclusively to indicate the magnetic field strength, while the subscript B is used 
to label the quantum numbers of monomer B. 

The plane waves can be expanded in terms of partial waves as 

00 L 

e^k-R^^^Y^ ^LJ^^kR)YLMQB.,<^R)YlUekAk), (10) 

L=0 M=-L 

where jiikR) is a spherical Bessel function of the first kind, the Yl^m functions are spherical 
harmonics, and the superscript * denotes complex conjugation. If we now substitute Eq. ([2]) 
for the particle interaction and use Eq. to describe the molecular asymptotic states, we 



7 



obtain 



La Li, 

Mt Ma q 



X 



[SaMHSbMH I [Sa ® 5^]^ |5^Ml"j, S^Mg), (11) 



where the last factor represents an integral over the spin coordinates. The integral over R 



can be performed analytically and gives, for ka < kb (see also Ref. 

lM^uAWR-l[^f 

r(^) 
r (i^i^^) r {La + 1) 

La- Lb-l La + Lb ^ 3 k^ 



X 2F, { + ^1, 3 , (12) 



where F is the Gamma function and 2F1 is Gauss's hypergeometric function. The integral 
over Qfi gives 



^L6,A4(^fl)C'2 -g(fiil)^L„,M„(^R)c?^^R 

^R 



V(2L. + 1)(2L. + 1)(-1)". I 2 2 L 

y -A4 A4 

with the terms in large round brackets denoting Wigner 3j symbols. The last 3j symbol 
readily implies that q = Ma — Mb. Finally, for the spin-dependent term we find 

{SaMI'ISbMI'I I [Sa ® Ss]^^^\SaMI'1 SsMfl) = 



V5{-1 



,q+SA+SB-M\"'-M 



3 A 



x[Sa{Sa + 1)(2^^ + 1)Sb{Sb + 1){2Sb + 1)] 



1/2 



^ I 1 1 2 \ / I Sa \l Sb 1 I 

^ ^ M'^M'^-q]\ -M£; M'^ Mfl ] \ -M12 M'^ M^] ' ' 



M'^M'b 



Note that the sums over and M'^ collapse for given values of M^^ and M^^, since the 
last two 3j symbols require that = M_^^^ — M^g^ and M'^ = M^^^ — M^^. Furthermore, 
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we have M^^ + = g so that M5 - M„ = M^"]^ + M^g] - M^^l - ^sl = ^^s- The sums 
over Mai and q in Eq. ( !TT|) therefore reduce to a single sum for any individual matrix 
element. The differential cross section is now readily calculated by substituting Eqs. ffTTl) - 
(p^ into Eq. (E]). 

The integral cross section is obtained by integrating daa^b/ dTlb over all orientations of 
the outgoing wave and averaging over all directions of the incoming collision flux, 

1 /" ,„ /■ d(Ja^b{^a) 



An 



dVta 



-dVth. 



(15) 



Using the orthogonality relation j j^i{yL)YL' M'iS^)d^ = Sl,l'Sm,m'' following 
expression for the BA cross section for {SAMg""^, SbM^"^) — )• {SaM^^, SbM^^) transitions 
induced by K„agn.dip: 



BA 



2%' 



La Li, 

xSAiSA + mSA + l)SBiSB + 1)(2^B + 1) 



1 



2La-l 



X 




^^i^^) r {La + 1) 

Lb — ^ La + Lb 



X 



X 



2 2 

Lf, 2 La 

-{Ma + AMs) ^Ms Ma 



AMs^ AMs^ 



' 2'k? 




Sa 1 Sa 
-M^l AMs, Mfl 




-1 2 



5., 
M, 



(a) 



(16) 



with AMs, = M^sl - M^l and AMs^ = Mf^ - M^^. The cross section for a specific 
incoming partial wave La and a certain outgoing wave L^ is obtained by simply omitting the 
sums over La and L^. We also point out that in the limit of ka ^ A;?,, which holds for ultracold 
exothermic collisions, the hypergeometric function 2F1 becomes 1 (see e.g. Eq. 15.1.1 of 
Ref. 0) 

and the energy dependence of the cross section is aa->,b ~ ^a^" ^ ~ E^°- 
The cross-section behaviour as a function of B is then, for B ^ E/ (gsfJ^sAMs), <Ja~^b ~ 
j^^-'^La ^ j^i/2-La_ jsjote that this i?-dependence is different from the threshold law derived 



by Volpi and Bohn 55|]. They considered the case of spin-changing transitions induced 
inside the centrifugal barrier of the exit channel, and found that the cross section behaves 
as aa-^b ~ 5^6+1/2. In our case, however, the spin- flip takes place at long range, outside 
the centrifugal barrier, and hence we find a different result. The long-range mechanism for 
dipolar spin relaxation will be addressed in detail in Sec. IIII A[ 

Equation (fT6!) is valid for any paramagnetic species that can be represented as in Eq. 
([9]). We note that, in the case of identical particles, the cross section must be multiplied by 
a factor of 2 if both monomers are in the same initial state, i.e. if Mg"^ = M^^^ (see e.g. 



Appendix B of Ref. [56|). This also applies to collisions between two magnetically trapped 
NH molecules, for which M^^^ = M^^l = 1. 



C. Analytical distorted-wave Born approximation 

As will be shown in Sec. Illlt the first-order BA is very accurate at low collision energies, 
but starts to deviate from the CC result at high energies and strong magnetic fields. One of 
the causes for this discrepancy is the phase shift in the incoming scattering channel. In order 
to quantify this effect, we have developed an analytical distorted- wave Born approximation 
in which the phase shift in the incident plane wave is explicitly included. 

Our starting point for the analytical DWBA is again Eq. (|6]), but now we approximate 
the incoming wave function jfei''''') as an elastically distorted wave \k'^), 

La Ma 

where h^^^ and hf^^ are spherical Hankel functions of the first and second kind, respectively, 
and Si^"^ is the elastic S-matrix element that contains the phase shift for the incident 
scattering channel {SaM^""^, SbM^^ , La, Ma) . The S*!^"^ matrix elements for NH-NH are 
taken from the full CC calculations described in Sec. Ill A[ The Hankel functions are defined 
in terms of regular and irregular spherical Bessel functions as 

hi\z) = jiiz) +iyL{z), 

hf\z) = jL{z)-tyL{z), (18) 

where yL{z) is a spherical Bessel function of the second kind. We note that the wave function 
of Eq. ([T7|) is unphysically divergent at the origin for nonzero phase shifts {Sia^^ ^ 1), but 
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matches the (exact) CC wave function at sufficiently large R. Hence, the approximation of 
Ifci"*"^) ~ constitutes an improvement over the ffist-order BA if the coupling occurs at 
long range. 

The radial part of Eq. f[T71) may also be written in terms of the transmission matrix 
element Tif^^ = i - si'^''\ 

\[hfsKR)+s^j^;^Sk^R)\ = \[hfskaR)+h'^^:^sk^R)-Tj^yh^;;^sk^R)^ 

= jLAKR)-lT^'J^^hil{KR). (19) 

Substitution of Eqs. f|T7j) and f|T9l) into ffTTl) for the matrix element over \4iagn.dip gives a 
radial integral of the form 

jL,ihR)^n^iKR)dR - irif") J jLAhR)^hilikaR)dR. (20) 



Note that the ffist integral is identical to that of Eq. (1121) . Using Eq. ( JTSl) . we may write the 
second integral of Eq. ( I20l) as 

jLAhR)^hiliKR)dR = 

J jL,ihR)^jLAkaR)dR + t J 3LAkbR)^yu{kaR)dR. (21) 

Again we observe that the ffist integral on the right-hand side is given by Eq. f lT2|) . The 
second integral on the right-hand side is convergent only for L;, > + 1 and gives, for 
ka < kh and integer La and L^, 

La + l 



J 3LAhR)^yLSKR)dR = (^f^ 
^ r (L. + 1) r (-^-^-^ 



-L„ - - 2 - L„ - 1 1 _ ]^ 

2 £ 



( , ^ ^ r , i - ^a, S ) . (22) 



We can now replace the radial integral in Eq. ( ITT]) by the expression of Eq. (120|) to obtain 
the matrix element over Knagn.dip in our distorted-wave Born approximation. Substitution 
into Eq. ([6]) gives the differential DWBA cross section, and Eq. (|T5l) subsequently yields the 
integral cross section. The final expression for the DWBA spin-inelastic cross section due 
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to K 



magn 



dip IS 



0" 



DWBA 



La Lf, 



2K 



xSa{Sa + 1){2Sa + 1)Sb{Sb + 1){2Sb + 1) 



X 



E 



Ma 



h 

La — Lb — 1 La + Lh 



La 



r (^^^) 



r( 



Ll, — La+3 



)r(i. + i) 



X 



r(L. + |)r( 



Lfj—La — ^ ' 



-Zvn — Lh — 2 Zvf, — Ln 



' 2 e 



Lb 2 La 

-(M„ + AMs) AMs Ma 



X 



X 




1 1 2 

AM5^ AM^^ AM5 



-\ 2 



5a 



5a 



5. 



-M'il AMs, Mfl 



-M'il AMs, Mfl 

(Ma) 



(23) 



The BA result of Eq. ( !T6|) is recovered in the hmit of Taa " — ?• 0. We emphasize that, in 
contrast to the BA, the sums over La and Lb in our DWBA expression should be restricted 
such that Lb > La + I [see Eq. ( 122|) ]. We also note again that, for indistinguishable particles 
such as NH + NH, the cross section must be multiplied by 2 if the monomers are in the 
same initial state. 



III. RESULTS AND DISCUSSION 



Numerical results 



We first discuss the numerical results for NH-NH obtained from full CC calculations. 



Previous theoretical work 39|] has shown that the intermolecular magnetic dipole interaction 
is the dominant trap-loss mechanism for NH-NH at low collision energies and small magnetic 
fields, while at higher energies and fields the intramolecular couplings become increasingly 
important. Here we will address only the intermolecular coupling term and provide a careful 
analysis of its contribution to the total inelastic cross section. 
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As explained in Ref. [39|, the contribution from Knagn.dip is most easily understood by 
considering the adiabatic potential curves. We will repeat part of this discussion here for the 
sake of clarity. Asymptotically, the adiabatic curves correspond to the molecular eigenstates 
Xi, and 7b, and at finite R each curve also contains a centrifugal barrier (see Fig. 1 of Ref. 



39|). Thus, at long range, the adiabats can be labeled by \^a,1b) and L, and therefore also 



correlate to scattering channels. It has already been noted in Refs. [3^ and j39| that several 
adiabatic curves are narrowly avoided due to the intermolecular magnetic dipole interaction, 
and that the spin-fiip induced by Knagn.dip takes place at the corresponding crossing. If we 
neglect the weak intramolecular spin-spin and spin-rotation couplings so that Eq. holds, 
we can define the avoided-crossing points Rc as 

2fiRi 

where La and denote the values of L for the adiabats correlating to the incoming and 
outgoing channels, respectively. The energies at which the crossings occur are given by 

= (25) 

defined relative to the threshold of the incident channel. We must point out that, since 
Knagn.dip coutaius a sccoud-rauk tensor in and first-rank tensors in the monomer spin 
coordinates, the avoided crossings occur only if La and Li, differ at most by 2 and M^g^^ and 
M^^ {i = A,B) each differ at most by 1. Thus, not all crossings are avoided. 

It can be deduced from Eq. fl24l) that, for small to moderate field strengths, the crossing 
points Rc are located at very long range. Therefore, the spin-change due to V^agn.dip can 
occur without having to overcome the centrifugal barrier in the outgoing channel. More 
specifically, the \Ms^ = Mgg = = 0) incident channel of NH-NH can couple with the 
\Msj^ = 1, Msg = 0,Lb = 2) and \Ms^ = Msg = 0,Lb = 2) outgoing channels even at zero 
collision energy. Hence, at low energies and relatively low magnetic fields, the intermolecular 
magnetic dipolar interaction is the main source of trap loss for NH-NH. 

It also follows from Eqs. (1241) and fl25|) that the curve crossings for higher partial waves, 
e.g. for the \Ms^ = Msg = = 2) — )■ \Ms^ = Msg = 0,Lb = 4) transition, become 

lower in energy as the magnetic field decreases. This implies that, at a fixed collision 
energy, the = 2 — )■ = 4 channel transitions open up below a certain 5-value. We 
will denote this critical magnetic field strength as B^. Figures [1] and [2] show the state- 
to-state inelastic NH-NH cross sections for different La — t- Lh channel transitions as a 
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function of B at collision energies of 10~^ K and 10~^ K, respectively. The Be values for 
the La = 2 — )■ Lfe = 4 and = 4 — Lf, = 6 crossings are also indicated. For the 
\Ms^ = = 1) — i- IM^^ = = 0) transitions, with AMs = 2, the numerical values 
are B^ = 8.67 x lO^^ G at 10"*^ K and B^ = 4.09 x lO^^ G at lO^^ K for L„ = 2 = 4, 

and Be = 8.67 G at 10"'^ K and B^ = 4.09 G at lO^^ K for L„ = 4 ^ L;, = 6. The critical 
field strengths for the \Ms^ = Ms^ = 1) ^ \Ms^ = l,Msg = 0) + \Ms^ = 0,Msg = 1) 
transitions, with AMs = 1? are twice as large as those for AMs = 2. It can be seen that 
the inelastic cross sections for = 2 — t- L5 = 4 and = 4 — > L^, = 6 indeed decrease 
as B exceeds the corresponding Be value. This S-dependence is remarkable, considering 
that higher partial waves typically contribute only if the exothermicity in the outgoing 
channel is large. Due to the long-range nature of the magnetic dipole interaction, however, 
the scattering of higher partial waves becomes increasingly important as the exothermicity 
decreases, because the crossings points Rc occur at a larger distance. 

The influence of the kinetic energy on the inelastic cross section can also be understood 
in terms of the adiabatic curve crossings. For a given magnetic field strength, the avoided 
crossings for = 2 — )■ L;, = 4 and = 4 — )■ = 6 are accessible only if the collision energy 
exceeds the value of Eq. f l25|) . It followed from Eqs. f l24|) and f l25|) that, if Ec increases, 
the critical field strength Be increases as well, and higher partial waves can contribute over 
an increasingly wide range of fields. This is also reflected in Figs. [1] and |2l In the ultracold 
regime, at a collision energy of 10"*^ K (Fig. [I]), the B^ values for = 2 ^ L?, = 4 and 
La = 4 —7- Lb = 6 are relatively small and the s-wave incident channel {La = — > L^ = 2) 
is strongly dominant at all field strengths above B ^ 10^^ G. At 10^'^ K, however, the Be 
values for the higher partial-wave channels are much larger, and we find that the L^ = 2 
and 4 incoming channels play a significant role at all magnetic field strengths below i? ^ 10 
G. A more detailed discussion on the energy dependence of the spin-inelastic cross section, 
based on the Born approximation, will be given in the next section. 

B. Comparison with BA and DWBA 

Before comparing our numerical results with the analytical BA and DWBA expressions, 
we must first point out that Eqs. ( 1161) and ( l23l) apply only to collisions in which Ms^ and Msg 
each change at most by 1 and L changes at most by 2. Furthermore, the integral of Eq. ( |22l) is 
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defined only if > La + 1, and tlie DWBA cross section of Eq. (l23i) is therefore valid only for 
La-> La + 2 transitions. Figure [3] shows the |M^"^ = 1, M^"-^ = 1, La) {Mf^, M^^.La + 2) 
cross sections as a function of 5 at a collision energy of 10~^ K. The cross sections are 
summed over all final states with \Mf^ - Mfl\ < 1 and \Mf^ - Mfl\ < 1. Figure H 
shows the results as a function of collision energy at a magnetic field strength of 1 G. It can 
be seen that the BA results are in very good agreement with the cross sections obtained 
from full CC calculations, in particular at low magnetic fields and low collision energies. 
At high fields and energies, the numerical cross sections exhibit several resonance features 
that arise mainly from the intramolecular spin-spin coupling term. Note that this coupling 
term is not included in the (DW)BA. Previous work has shown that the intramolecular spin- 
spin coupling becomes increasingly important as the kinetic energy in the outgoing channel 
increases, and, for B > 10^ G and E > 10~^ K, causes almost the same amount of spin 
relaxation as the intermolecular magnetic dipolar interaction 



. Hence, the BA result of 



Eq. ( 1161) deviates from the full CC result at high energies and field strengths. 

It can also be seen in Figs. |3] and H] that the analytical distorted-wave BA cross section, 
which contains an extra term due to the phase shift in the incoming channel, is in slightly 
better agreement with the numerical = — > Lf, = 2 cross section than the BA result. In 
particular. Fig. |3] shows that the BA cross section for = — > L;, = 2 starts to deviate 
from the CC calculations around B ^ 1 G, while the DWBA is accurate up to i? ^ 100 G. 
Thus, in the region between 1 and 100 G, the inelastic Lq = — ?■ Lf, = 2 cross section can 
be completely attributed to the intermolecular magnetic dipole interaction and to the phase 
shift in the incident channel. For the higher partial-wave channels, however, the analytical 
DWBA cross section deviates significantly from the full CC result at high fields and low 
energies. This is due to the (/cfe/fca)^""*"^ term in the expression for a^^^^ [Eq. f l23|) ]. which 
tends to infinity if ka <^ kb- Even for very small phase shifts, this term will dominate the 
DWBA inelastic cross section for > if the collision energy is small and the exothermicity 
is large. More specifically, we estimate from Eq. ( I23l) that the DWBA cross section diverges 
if {ka/kby^"-^^ ^ Ta^''\ and hence the effect is most pronounced for large La- We point 
out that the origin of the (kb/ka)^"-^^ term lies in the irregular spherical Bessel function 
ULai^aR) [Eq. ( 122|) ]. which enters the asymptotic wave function if the phase shift is nonzero. 
At short range, the yi^ function tends to infinity and ultimately leads to the unphysical 
behaviour observed in Figs. E] and HI A possible remedy for this problem is to evaluate the 
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integral of Eq. (l22l) only for i?- values larger than a certain cutoff radius. However, such an 
approach requires careful numerical analysis and falls out of the scope of the present study. 
Nevetheless, based on the results shown in Figs. [3] and HI we conclude that the BA gives more 
meaningful results than the analytical DWBA at most of the energies and fields considered 
in this work. As a final point, we note that the numerical phase shifts for the higher partial- 
wave channels are orders of magnitude smaller than the s-wave scattering phase shift, and 
our DWBA results would not be substantially improved by including a phase shift in the 
outgoing channel. 

As derived in Sec. II1B[ the threshold behaviour of the BA spin-inelastic cross section 
in the limit of ka ^ is given by aa^b ~ b^/'^^l^ aa^b ~ E^'^^^^'^. Indeed, we find 
that the = — j- L;, = 2, = 2 — )■ = 4, and = 4 — )■ Lf, = 6 inelastic cross 
sections at 10~^ K behave as i?^/^, B'^^"^, and respectively, for field strengths above 

B ^ 5 X 10^^ G (see Fig. |3]). Similarly, the cross sections at 5 = 1 G follow an E'^^^, E^/^, 
and E'^/'^ dependence, respectively, at collision energies below E ^ 10~^ K (see Fig. Hj). The 
validity regions of these threshold laws can also be explained in terms of the Knagn.dip-induced 
avoided crossings discussed in Sec. 1111 A[ The critical magnetic field strengths below which 
the = 2 — 7- Lj, = 4 and = 4 — Lf, = 6 crossings are energetically accessible are on 
the order of Be ~ 10^^ G for a collision energy of 10~^ K (see Fig. [1]). If the magnetic field 
strength exceeds B^, the crossings for the higher partial-wave channels are inaccessible and 
the corresponding scattering process can proceed only by (non-classical) tunneling through 
the centrifugal barrier. Hence we find the quantum-mechanical threshold behaviour at fields 
above B ^ 10^^ G. For field strengths below Be, the approximation of ka <^ kb breaks 
down and the 5-dependence follows from the explicit evaluation of Eq. (fT6|) . That is, the 
5-dependent threshold behaviour for B < Be is given by {ka/kbY^"'^ multiplied by the 
hypergeometric function. If the magnetic field is so small that ka ~ fcb, the field dependence 
becomes negligible and the cross section flattens off to a constant value. In order to explain 
the energy dependence, we apply Eqs. (1241) and (1251) to determine the lowest possible Ec 
values at which the avoided curve crossings can occur. At a magnetic field of 1 G, the 
corresponding values are Ef. = 5.8 x 10^^ K for the = 2 — t- Lf, = 4 transition and 
Ec = 1.2 X 10~^ K for = 4 — )■ L;, = 6. Since the crossings for the higher partial- wave 
channels are inaccessible if E < E^., we recover the quantum-mechanical threshold law at 
collision energies below E ^ 10~^ K. 
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The results presented so far apply only to collisions where M^^ and Msg decrease at 
most by 1 and L increases by 2. The total spin-inelastic cross section, however, contains 
contributions from all (symmetry-allowed) outgoing partial waves and all final states, i.e. 
also the states with M^'^ = —1 and = —1. Let us now compare the BA and DWBA 
results with the numerical total spin-inelastic cross sections for magnetically trapped NH 
{M^sl = ^£ = 1)' summed over all possible incoming partial waves and all outgoing 
channels. The total BA cross section is obtained by performing the sums over La and L;, in 
Eq. ( !T6|) for all possible (even) partial waves. To calculate the total DWBA cross section, we 
perform the sums over La and L;, in Eq. f l2^ for all possible (even) La values and L^ = La + 2. 
Since the numerical scattering calculations were carried out for Lmax = 6, we also took this 
maximum value for La and Lf, in the (DW)BA expressions. 

The total inelastic cross sections are presented in Figs. [5] and [61 It can be seen that the 
BA is generally in much better agreement with the full CC result than the DWBA, except 
for a small region near 10 G at 10^^ K (Fig. E]) and near 10^'^ K at 100 G (Fig. [6]), where 
ka ~ kb- As noted previously, the deviation of the DWBA at high B and low E is due to 
the (kb/ka)^"'^^ term in Eq. ( I23i) . which causes unphysical behaviour if kb ^ ka- At low 
magnetic fields and relatively high energies, in particular at E = 10^'^ K (see Fig. [5]), we 
find that the DWBA cross section also deviates from the numerical result. In this region, 
the spin relaxation arises mainly from the La = 2 Lb = 2 transition, and, to a smaller 
extent, also from the La = Lb = 2 and La = 2 ^ Lb = transitions. Since the 
total DWBA cross section is restricted such that Lb = La + 2, the most dominant inelastic 
transitions at low B are thus not included in the DWBA. Similarly, the total DWBA cross 
section as a function of energy (Fig. [6]) shows a clear discrepancy with the numerical result 
at B = 0.01 G for nearly all energies considered, and at 5 = 1 G for E > 10~^ K. This is 
also due primarily to the missing La = 2 ^ Lb = 2 channel transition. 

It can also be seen that the total BA cross section, which does include all possible La — )■ Lb 
transitions, agrees over a much wider range of E and B, but deviates from the full CC result 
at high fields and high collision energies. As already discussed in the first paragraph of this 
section, the deviation partly arises from the intramolecular spin-spin coupling term, which 
contributes significantly to the numerical cross section as the kinetic energy in the outgoing 
channel becomes large. Moreover, the BA only includes contributions from final states with 
Ms^ = 0, 1 and Msg = 0, 1, while the total numerical cross section also contains terms 
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with Msj^ = — 1 and Msg = —1. Previous work has shown that, as the intramolecular 
spin-spin coupling term becomes increasingly irnportant, the state-to-state cross sections for 



Uf^ = -1 and/or M^'l = -1 increase as well [39|. Although the dominant mechanism for 



M^g"^ = 1—7' M^^ = — 1 transitions is likely to be the intramolecular spin-spin term, which 
can decrease Ms- by 2 directly in first order, the intermolecular magnetic dipolar coupling 
term may also induce such spin-changing collisions in second order. This effect is included 
only in the full CC calculation, and hence this may represent another source of discrepancy 
between the BA and the numerical result. 



IV. CONCLUSIONS 

We have presented a detailed theoretical study on the role of the intermolecular mag- 
netic dipole-dipole interaction in (ultra)cold collisions of magnetically trapped NH(X^S~) 
molecules. The inelastic cross sections for Zeeman relaxation have been obtained from rigor- 
ous coupled-channel calculations and from analytical results based on the (distorted-wave) 
Born approximation. The derived expressions for the analytical cross sections are valid for 
any two paramagnetic species for which the electronic spin and its space-fixed projection 
are (approximately) good quantum numbers, but we have applied them only to the NH -|- 
NH system. 

We have found that the scattering of different partial waves, induced by the magnetic 
dipolar coupling, is most easily understood by considering the adiabatic potential curves. 
The intermolecular dipolar coupling term induces avoided crossings between certain adiabats 
at long range, which in turn may lead to Zeeman relaxation. The cross-section behaviour as a 
function of energy and magnetic field is, to a large extent, determined by whether or not these 
avoided crossings are energetically accessible. Remarkably, the avoided crossings for higher 
partial waves become lower in energy as the magnetic field strength decreases, implying 
that the corresponding channels open up below a certain critical field strength. Indeed, it 
was found that the scattering of higher partial waves becomes increasingly important as the 
exothermicity decreases. 

The validity regions of the analytical BA and DWBA have been determined by comparison 
with numerical close-coupling calculations. We have found that the BA is accurate over a 
wide range of collision energies and fields, but starts to deviate from the numerical cross 
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sections at energies above 10~^ K and fields above 10^ G. This is mainly due to 
the effect of the intramolecular spin-spin coupling term, which is neglected in the BA but 
contributes significantly to the numerical cross section as the kinetic energy in the outgoing 
channel becomes large. The analytical distorted-wave Born approximation, which accounts 
for a phase shift in the incident channel and thus represents a correction to the BA, gives 
more accurate results than the BA in the case of s-wave scattering. For higher partial-wave 
scattering, however, and in particular at high magnetic fields and low energies, the DWBA 
cross section shows unphysical behaviour and diverges to infinity. Furthermore, the derived 
DWBA expression is valid only for collisions where the partial-wave angular momentum 
is increased by 2, while the total numerical cross section contains contributions from all 
possible outgoing partial waves. More specifically, at fields below 1 G and energies near 
10~^ K, the dominant contribution to the inelastic cross section is the La — 2 ^ Lj, — 2 
channel transition, which is not included in the DWBA. Hence we conclude that the BA, 
which contains all possible partial-wave contributions and does not show any unphysical 
behaviour, is generally in much better agreement with the full GG result than the DWBA. 

Although we have focused only on NH('^S^) + NH(^E~) collisions in this study, the theory 
and main conclusions should be general for any two (ultra) cold paramagnetic species. 
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FIG. 1: State-to-state collision cross sections for the \Msa = ^^^Sb = 1) ~^ \^Sa — ^-i^Sb — 0) 
and \Ms^ = l,Msg = 1) ^ IM^^ = l,Msg = 0) + \Msj, = 0,Msg = 1) transitions of NH + 
NH, obtained from full CC calculations as a function of magnetic field at = 10^^ K. The (blue) 
lines marked with squares correspond to = — t- L;, = 2 transitions, the (red) lines marked with 
triangles correspond to La = 2 — t- L(, = 4, and the (green) lines marked with circles correspond 
to -Lq = 4 — )• Lfe = 6. The vertical lines indicate the Be values below which the crossings for 
La = 2 — )■ Lf, = 4 and = 4 — t- Lj, = 6 are energetically accessible. 
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FIG. 2: State-to-state collision cross sections for the \Msa = ^^^Sb = 1) ~^ \^Sa — ^-i^Sb — 0) 
and \Ms^ = l,Msg = 1) ^ IM^^ = l,Msg = 0) + \Msj, = 0,Msg = 1) transitions of NH + 
NH, obtained from full CC calculations as a function of mag netic field at ^ = lO"^ k. The (blue) 
lines marked with squares correspond to = — t- L;, = 2 transitions, the (red) lines marked with 
triangles correspond to La = 2 — t- Lf, = 4, and the (green) lines marked with circles correspond 
to Lq = 4 — )• Lfe = 6. The vertical lines indicate the Be values below which the crossings for 
La = 2 — )• Lf, = 4 and = 4 — t- Lj, = 6 are energetically accessible. 
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FIG. 3: Total spin-inelastic collision cross sections for the jM^^ = l,Msg = !)—)• \Ms^ = 
0,Ms^ = 0) and \Ms^ = l,Ms^ = 1) ^ \Ms^ = l,Ms^ = 0) + \Ms^ = 0,Mss = 1) transitions 
of NH + NH, calculated as a function of magnetic field at 10^^ K. Different colors correspond to 
different La — s- Lf, channel transitions. 
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FIG. 4: Total spin-inelastic collision cross sections for the jM^^ = l,Msg = !)—)• \Ms^ = 
0, Ms^ = 0) and \Ms^ = 1, Ms^ = 1) ^ \Ms^ = 1, Ms^ = 0) + \Ms^ = 0, Ms^ = 1) transitions of 
NH + NH, calculated as a function of collision energy at a magnetic field strength of 1 G. Different 
colors correspond to different La — s- Lf, channel transitions. 
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FIG. 5: Total spin-inelastic collision cross sections for two magnetically trapped NH molecules, 
calculated as a function of magnetic field. Different colors correspond to different collision energies. 
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FIG. 6: Total spin-inelastic collision cross sections for two magnetically trapped NH molecules, 
calculated as a function of collision energy. Different colors correspond to different magnetic field 
strengths. 
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